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Abstract 

This paper considers the issue of Bose-Einstein condensation in a weakly 
interacting Bose gas with a fixed total number of particles. We use an old 
current algebra formulation of non-relativistic many body systems due to 
Dashen and Sharp to show that, at sufficiently low temperatures, a gas of 
weakly interacting Bosons displays Off-diagonal Long Range Order in the 
sense introduced by Penrose and Onsager. Even though this formulation is 
somewhat cumbersome it may demystify many of the standard results in the 
field for those uncomfortable with the conventional broken symmetry based 
approaches. All the physics presented here is well understood but as far as 
we know this perspective, although dating from the 60's and 70's, has not 
appeared in the literature. We have attempted to make the presentation as 
self-contained as possible in the hope that it will be accessible to the many 
students interested in the field. 
Pacs numbers: 03.75.Fi 05.30.Jp 
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I. INTRODUCTION 



The upsurge of activity in recent years in the area of Bose Einstein Condensation (BEC) 
in dilute Bose gases has led to remarkable experimental advances and innovations unimag- 
inable ten years ago when the first condensate in trapped alkali atoms was reported [|IJ. In 
parallel with these beautiful experiments a better theoretical understanding is also emerg- 
ing: These systems are the first experimental realization of the weakly interacting Bose 
condensed gas studied long ago by Bogoliubov, Beliaev, Lee and Yang, and Gross and 
Pitaevskii 0. Moreover, there are issues specific to the trapped alkali atoms - the multi- 
component structure of the condensate due to the internal atomic hyperfine structure, the 
strongly inhomogeneous nature of the condensate reflecting the strong inhomogeneity intro- 
duced by the trapping potential, the rather small number of particles involved in some cases 
(as low as 10 5 compared to the typically 10 23 in liquid helium), as well as the metastability 
and sometimes instability of the gaseous phase with respect to self-binding into liquid or 
solid droplets - which make these systems rich and interesting in their own right. In fact, it 
is probably fair to say that a first principles theory, especially as far as dynamical properties 
are concerned, is not yet fully developed and tested 0. 

This paper was motivated by a heated discussion between the atomic physics and quan- 
tum optics communities on one hand and the condensed matter physics community on the 
other, concerning the necessity of using the notions of broken symmetry and order pa- 
rameters in describing the physics of BEC in trapped gases. Indeed, with the exception of 
variational ground state wavefunction approaches and computer simulations, all quantitative 
many-body treatments of BEC rely on the notion of a broken symmetry in which the ground 
state of the system breaks a global symmetry of the original Hamiltonian In the partic- 
ular case of BEC this is the global U(l) gauge rotation generated by U(0) = exp (^—iON^j 
associated with the conservation of the number of particles, [H, N] = 0. Formally this is 
implemented in explicit calculations by using the so-called "symmetry broken ensemble" 0] 
characterized by a density matrix, p, with the property, [p, N] ^ 0, in spite of particle 
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number conservation by the Hamiltonian. Within this description the BEC displays a finite 
"order parameter", /^(r)^ 7^ (/^d(r)\ 7^ 0), where ^(r) (\&t(r)) is the Boson annihilation 
(creation) operator at position r, and (..) denotes the ground state expectation value (at zero 
temperature, T = 0) or the expectation value in the appropriate thermodynamic ensemble 
(for T^O). One should imagine that our Bose gas can extract particles from or inject par- 
ticles into a sufficiently large particle reservoir disturbing the strict number conservation; 
the Bose condensate would then correspond to a coherent state built from of superposition 
of states with different numbers of particles allowed by the contact with the reservoir. 

It is no surprise that this formal and unintuitive construction appeared somewhat myste- 
rious to many unfamiliar with superfluid liquid Helium physics who felt uneasy in applying 
these ideas to an isolated system with a fixed number of particles such as the trapped 
gases. This led to a number of reformulations and generalizations of the Gross-Pitaevskii 
and Bogoliubov-de Gennes equations which respect the U(l) symmetry ||. These schemes 
are more awkward than the conventional many-body techniques, less amenable to analyti- 
cal approximations and more difficult to extend to finite temperatures and systems far from 
equilibrium. Here we introduce the reader to yet another [/(l)-symmetric approach based on 
the current algebra formulation of non-relativistic many-body systems initiated in the 60 's 
by Dashen and Sharp || . As all other particle conserving schemes this also has its unsettling 
features (like, as shown below, the appearance of singular operators); nevertheless, as the 
approach only involves operators which commute with the total number of particles, it can 
be studied with well-known field theory techniques and can be easily extended to treat finite 
temperatures and non-equilibrium situations. We note that the current algebra descrip- 
tion can be interpreted as the operator version of the [/(l)-symmetric functional integral 
representation due to Popov J/J. Different aspects of Bose condensation using the Dashen- 
Sharp current algebra formulation have been discussed in the past || but, as far as we know 
ours is (i) the first discussion of inhomogeneous systems and (ii) the first computation of 
single-particle correlation functions. 

Historically, the description of BEC in systems with a fixed number of particles follows 
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from the realization || that, in the thermodynamic limit (with iV — > oo with the volume 
Q —>■ oo with n = N/Q fixed), the Bose condensed state displays long-ranged correlations 
in the single particle density matrix, p(r; r') = ^"'"(r)^'(r / )^>, where ^(^(^(r)) is the 
creation (annihilation) operator for a boson at site r = (x,y,z). More precisely, as p(r, r') 
is hermitian with respect to the r and r' indices, it can be expanded in terms of its complex 
eigenfunctions, <f> a (r) and real eigenvalues, A a as, 

p(r;r') = (^(r)^')) = £ A a C(r)0 Q (r'). (1) 

a 

BEC in the noninteracting gas is signaled by a macroscopic eigenvalue associated with the 
zero momentum state, Ak =0 = Uq (T), where in this case a labels momentum eigenstates 
0k(r) = [exp(ik • r)]/VTi, n\ (T) is the condensate density of the noninteracting gas at 
temperature T and Q is the volume of the system. BEC displays long range correlations 
in the sense that lim\ r ^ r i\^ 00 p{r\Y l ) = n[ ) O ' ) (T)^ =o (r)0k=o( r ')5 while the contributions from 
higher momentum states oscillate away in the limit. By adiabatic continuity, condensation 
in the interacting Bose gas is defined by the presence of one eigenfunction of the single 
particle density matrix, 0o( r ) [11 |> with a macroscopic eigenvalue, Ao = uq. It's not hard 
to show that, at T = 0, $(r) = y/no4>o(r), often referred to as "the macroscopic wave- 
function" satisfies a non-linear Schrodinger equation, the time independent version of the 
Gross-Pitaevskii equation 0. In this case the system is said to display Off-diagonal Long 
Range Order (ODLRO) but (^(r)) = and thus, strictly speaking, the symmetry remains 
unbroken |fID[] 

Below we use the current algebra approach in a Bose gas at zero temperature to calculate 
the low energy excitation spectrum and to check for the presence of BEC in the ground state, 
in the sense of Penrose and Onsager B . 



II. THE CURRENT ALGEBRA APPROACH 



The discussion that follows is based in its entirety on the work of Dashen and Sharp || . 
The idea is to represent the Hamiltonian and all operators of the theory in terms of the den- 
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sity, n(r) = \I/t(r)\I/(r) and current g(r) 
which obey the current algebra: 



i/2M)[^ t (r)V^(r) - W f (r)tf(r)] operators 



[n(r), n(r')] = 



n(r),g(r 



M 

% 

M 



V r (5(r - r')n(r)) 
V^(r-r')g Q (r)) 



V?,(5(r-r')gV)) 



(2) 
(3) 
(4) 



(Hereafter = x,y,z index vector components, U — 1, and the hat will differentiate 
operators from classical fields whenever ambiguities can arise.) 

The main step of this formulation is to rewrite the Bose Hamiltonian in terms of n(r) 
and g(r) by using unity in the form, I = \&(r)[l/n(r)]\E'"'"(r) (H — 1) [12 

1 



H 



f 1 



W f (r 



n(r) 



V*(r)) + ^Mu; |r| 2 n(r) + |n 2 (r) 1 dr (5) 



8M 



Vn(r) - 2zMg(r) 



1 



n r 



Vn(r) + 2iMg(r)l + -Ma;o|r| 2 ri(r) + ^n 2 (r) I dr, 



where for simplicity the interparticle potential was replaced by the s-wave pseudopotential, 
v(r—r') = gS(r—r / ), with the strength, g = Aira/M, written in terms of the s-wave scattering 
length, a; and use was made of the identities, (V^(r)$(r)) = (Vn(r) — 2zMg(r))/2 and 
(^(r)V^(r)) = (Vn(r) + 21Mg(r))/2. We have included an external harmonic potential 
(with frequency Uq) which traps the particles in a finite region of space. 

Representations of the current algebra @, (||) and @ have been discussed extensively 
in the 70's It was shown that in an irreducible iV-particle representation of the current 
algebra matrix elements of operators such as (Vn(r) — 2zMg(r)) and (Vn(r) +2iMg(r)) are 
proportional to n(r), and thus the singular operator, l/n(r), in (f|), disappears in physical 
matrix elements. Below we will leave aside all rigor and manipulate expressions involv- 
ing l/n(r) formally with the expectation that, even if intermediate states of some of the 
calculations are ill-defined, the final answer is physically meaningful. 
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A. The Excitation Spectrum 



To obtain the low lying excitation spectrum we use a mean field approximation and 
expand the Hamiltonian (0) to quadratic order in fluctuations of the density and current 
around their values in the ground state. Below we only describe the more general case of 
the nonuniform system and extract the homogeneous gas results as a special limit. 

In the presence of the harmonic potential the ground state, |0), is characterized by 
an inhomogeneous particle density, (Q\n(r)\Q) = nc(r), and a (particle) current density, 
0|g(r)|f2) = gG( r ) ; to be determined by minimizing the mean-field energy functional, 



£csK?(r), <? G (r)] 



[(Vn G (r)) 2 + 4M 2 (g G (r)) : 



+ -Mcu \r\ 2 n G (v) + ^-n 2 G (r) } dr, (6) 



8Mn G (r) 2 "' 2 

subject to the fixed particle number constraint, / drn G (r) = N. As usual this con- 
straint is enforced by adding a chemical potential like term to @, Ecs[n>G(. r ), 9G( r )) — * 
Ecs[n G (r) , g G (r)} — \i j <irn G (r). The resulting ground state carries no current, g G (r) — 0, 
and a nonuniform density satisfying, 

(V r n G (r)) 2 



8Mn G (r) 







(7) 



The low lying, large length scale excitations of the system are described by the effective 
Hamiltonian, Hx, obtained from by making the replacement n(r) = n G (r) + fj(r) and 
keeping the terms leading (i.e., quadratic) order in the excitation operators, fj(r) and g(r): 



H 



x 




V r r)(r)) 2 + 4M 2 (g(r)) 1 



+ 



V r n G (r) 



8Mn G (r) \ // G (r; 
fj(r) (V r r)(r) + 2zMg(r)) + (V r f)(r) - 2zMg(r)) fj(r) 



+ § 



f) 2 {r) 



8Mn 2 G (r) 

For the purpose of obtaining the linear excitations, the current algebra (H), @ and (|) is 
replaced by the (linearized) approximation, 



[n(r),n(r')] = 0; 
f)(r),g(r' 



^V r (<5(r-r> G (r)). 



(9) 
(10) 



With these simplifications the equations of motion for the excitation operators, read: 



0_ 

Of 

d_ 

Ot 



f)(r,t) 
g(r,t) 



- V r -g(r,t) 



AM 2 



f)(r,t) 



Vi(l/n G (v))-A V r l/Jn G (r 



(11) 



(12) 



where summation over repeated indices is implied. Naturally the density fluctuations satisfy 
a continuity equation (|TT|) . 



For illustration consider the Thomas-Fermi (TF) limit [14J] in which one ignores the 
gradient terms in (0) and in the curly bracket in (|12|) . In that case, gn G (r) ~ Mu 2 (R 2 — 
|r| 2 )/2, and Equations ( |IT|) and (|T2|) combine to give, 

d 2 1 



(13) 



(R is the spatial extent of the condensate which, for a spherical trap with N particles, is 
given by R = a^o (lBNa/aho) 1 ' 5 in terms of the harmonic oscillator length, a^o and the 



s-wave scattering length, a.) This is precisely the equation first discussed by Stringari fl5 
and extensively studied since for a variety of trap geometries |J. 

It is clear that to obtain the results of the uniform system one must go beyond the TF 
approximation. This can be done either directly from equations (0), (|TT|) and (0) or by 
making contact with the conventional non-conserving Bogoliubov approach fl6|| . We will 
take the latter route. 



1. Correspondence to the Bogoliubov Approach 

We proceed by first noticing that the identity, V r (<5(r — r')n G (r)) = — n G {r') V r '5(r — r') 
suggests rewriting the current operator as g(r) = n G (r)V r <p(r) / 'M , where the operator 
0(r) is the canonically conjugate momentum density to r/(r): [f)(r),ip(r')] = i5(r — r') and 
[0(r) , (p(r')\ = 0. As in the example of the simple harmonic oscillator we can then make 
linear combinations of "coordinates" (17 (r)) and "momenta" (y3(r)) to construct creation 
(fet(r)) and annihilation (6(r)) operators. Here we choose, 
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n G (r)6(r) = - (fj(r) + 2in G (r)0(r)) 
S(r),S t (r')l =5{r-r'); [S(r),S(r')l = \&{r), &(r') 



0. 



It is not hard to see that, in terms of the new canonical Bose fields (H) becomes, 



H 



x 



dr 



V r 6 f (r)) • (V r 6(r 



n G (r) (S(r) + & f (r 



\S)i f (r) = g (*K r ) - 2m G (r)y3(r)) (14) 

(15) 
(16) 

_2MV'--V V-v-V ' 2 — v-/ v-v-/ ■ - v/y j p ( 17 ) 

where the only term in ([§D incorrectly reproduced in ([17|) are those involving the V r nG(r) 
terms proportional to r) 2 (r). To recover these terms and to eliminate unwanted terms pro- 
portional to 2 (r) requires going beyond the linear approximations ( |i"4|) and ([L5|) in the 
representation of 6(r) and 6^(r) in terms of 77(1") and <£>(r). (Note that all the problematic 
higher order terms involve gradients of n G (r) and thus vanish for a uniform system.) 

The Hamiltonian ([17]) is identical to that derived by Bogoliubov in the presence of a 
finite condensate order parameter ||16|| ; for the case of a fixed number of particles the order 
parameter vanishes and Hx should be interpreted not as a quasi-particle Hamiltonian but 
as the Hamiltonian describing the low lying density and current excitations of the system at 
a fixed total particle number. It is then no surprise that (|l7|) does not conserve the number 
of bosons, / drff(r)b(r). 

The quadratic form Hx is easily diagonalized by the Bogoliubov transformation, b(r) = 
En (u n (r)$ n (r) - <(r)/?t( r )) and &t(r) = En (<(r)/t( r ) " v n (r)$ n (r)), with the functions 
u n (r) and v n (r) satisfying, 



2M 
V 2 
2M 



+ gn G {r)J u n (r) 
+ gn G (r))v n (r) 



gn G (r)v n {r) 
gn G (r)u n (r) 



E n u„(r) 



- E n v n (r) 



(18) 
(19) 



and the orthonormality condition, / dr (w* (r)w n /(r) — u*(r)t; n /(r)) = 8 n>n r. 

For a uniform system, « k (r) = w k e lk ' r , f k (r) = f k e lk ' r , n G = n = N/Q, and 

these equations immediately lead to the well known Bogoliubov results |16[], u k = 

r 1 ~ r l 1 1 

^ \E^ 1 (e k + gn) + lj 2 and v k = \E^ 1 (e k + gn) - lj 2 , where £ k = ye k (e k + 2gn) 

(e k = |k| 2 /2M) is the Bogoliubov quasiparticle energy. 



B. Off-Diagonal Long Range Order 



To check for the occurrence of BEC in the sense of ODLRO is more involved as we 
need to calculate the single particle density matrix, p(r + ix, r) = ^A(r + ix;r)^) where 
A(r + ix; r) = \l/t(r + ix)\E'(r). At first sight it is hard to see how one might do this 
computation in an approach in which the basic variables are number and current densities. 
The solution can be found in the old work of Grodnik and Sharp ||. 

Without loss of generality we take the separation vector between the two points, r and 
r', in p(r,r'), along the x-axis. We then proceed by considering the differential equation, 



V x A(r + ix; r) = V x ^(r + ix)^(r) = V^^r + ix)^(r + ix) 



1 



_ ijrt( r _|_ ixW( r ) 

n(r + ix) 



- [V x h(r + ix) - 2iMg B (r + ix)] - . A(r + ix; r) 
2 n(r + ix) 



(20) 



where in the last two steps we used the resolution of the identity, I = ty(r)[l/n(r)]W(r). 
Together with the initial condition, A(r+ix; r)| x=0 = n(r), and the linearized approximation, 
[V«/n(r + ix') - 2tMg x {r + ix')] ^-tj ~ (St(r + ix')/ y/n G {r + ix')) , Eq. © leads 
to a "x-ordered" exponential which allows us to write down the general formula for the 
spatial correlation function: 



p(r + ix, r) = p(\x\) = (T x > \ exp [ dx'\ [V x >n(r + ix') - 2iMgJr + ix')] — — 

\ { Jo 2 n(r + ix') 



n r 



n G (r + ix)n G {r) 



exp 



St(i 



ix 



St(i 



S(r) +St(r) 



'n G [r + ix) yn G {r) 

The explicit computation of fl2~ID is then carried out by transforming to Bogoliubov 
quasiparticles and using the "disentangling" Baker-Hausdorff formula, exp(A + B) = exp — 
h[A, B]expAexpB (where [[A, B],A] = [[A,B],B] = 0) together with the condition that 
/3(r) annihilates the ground state, (3(r)\Q) = 0. The resulting expression, now written for 
an arbitrary separation vector, R, is expressed in terms of the Bogoliubov amplitudes, 1 ^ (r) 
and v n (r), as: 



(21) 
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I>(\H\) = jexp 

i-EP 



1 <(r + R) <(r) 

2.1. U G (r + R) Jn G (r) 



v„(r + R) 



x 



n G (r + R) \/rc G (r 



n G (r + R)n G (r). (22) 



n G (r) J \^n G (r + R) \Jn G (r) 
A similar calculation can be done for time dependent single particle correlation func- 
tions. For example, the local, time- dependent single-particle density matrix, p(r, t; r, 0) = 
A(r, t; r, 0)), where A(r, t; r, 0) = ^ f (r, t)^(r), is given by: 

1 1 



p(r, t; r, 0) = (T T j exp ^ dr^ <9 r ri(r, r) + (d T &(r, r)$(r, r) - ^(r, r)d T ^f(r, r 



nir 



n r, t 



n G (r)e^(T T 



* (ld r fj(r,T) 
exp I dr \ — — to T (p(r } r) 



n G (r e 



ifit 



exp 



o \2 n G (r 
St(r, t)-St(r) 



1 + 



nc(r) 



(23) 



6(r) +St(r ) 



The last two equalities in ( f23|) involve the linear approximation defined by ([14]) and (|15|). 
As with (|22|), the final answer is expressible in terms of Bogoliubov amplitudes, u n (r) and 
v n (r) as: 



p(r, t; r, 0) = n G {r)e %li I exp 



n G ir 



x 



x 



j + v- ^n(r) «(r) -<(r)) ^ _ 



(24) 



It is easy to see that in the homogeneous case (22) and (p3| ) become: 



^|RHooP(|R-|) =e %lxt lim t -+ 00 p(r, t; r, 0) 

2fi k i? k/ 



n 



2fin 



(25) 

k " ^ 

A number of features of ( p5"D are worth noting: (i) as expected the condensate density is 
equal to the total density in the noninteracting limit; (ii) expanding to leading order in g 
gives the same result for the depletion of the condensate as calculated from the conventional 
Bogoliubov approach |I6] (note that the ultraviolet cutoff required due to the linear approx- 
imation cancels to leading nontrivial order in g); (iii) in a system with a fixed number of 
particles the phase of the condensate contribution precesses uniformly at a rate determined 
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by the chemical potential, /i, in agreement with the Josephson relations ||17|| ; (iv) finally, 
as expected on general grounds a T = condensate occurs in two and three dimensions, 
whereas in one-dimension the single-particle correlation function decays algebraically both 
in space and in time (due to the infrared logarithmic divergence in the exponent). An ex- 
plicit analysis of the nonuniform Bose gas, including a discussion of the meaning of ODLRO 



in finite systems and a comparison with the results of Ma and Ho fllSfl , is left to a future 
publication. 

III. CONCLUSIONS 

We have considered the Bose gas in a general inhomogeneous potential within the current 
algebra approach to non-relativistic many-body systems due to Dashen and Sharp ||. Not 
surprisingly, we arrive at the same physics as that of broken symmetry approaches. The 
differences are somewhat subtle: strictly speaking, in the broken symmetry case the order 
parameter field displays collapse and revivals due to fluctuations in the particle number || . 
As a result, the Gross-Pitaevskii and Bogoliubov-de Gennes equations are adequate only for 
times short compared to the collapse time scale which becomes infinite only in the thermo- 
dynamic limit (in the uniform Bose gas r co u apse ~ y/~N). This effect was understood already 
in the 50 's through Anderson's classic discussion of broken symmetry in quantum antifer- 



romagnets fll9 |, for which, as in Bose condensates, the order parameter is not a constant of 
the motion. In all such systems, apart from the Goldstone modes present as a result of the 
broken symmetry, there exist modes - the "phase diffusion" mode in the case of BEC - with 
a frequency which vanishes in the infinite volume limit faster than that of the lowest Gold- 
stone mode. Our fixed-iV formulation only includes the physics of a fixed N sector within 
which the phase precesses at a constant rate given by the chemical potential. Collapse and 
revivals can be obtained only by averaging over systems with different values of N, as would 
be the case if our condensate was brought into contact with an ideal particle reservoir. 
Controlled {/(l)-symmetric techniques are especially important in treating small systems 
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and systems far from equilibrium where spurious dynamics of the condensate in symmetry 
broken ensembles may confuse some of the important physics. We expect that the theoreti- 
cal frameworks presented here and in |5| can be used to analyze the feasibility of a number 
of novel experiments on phase coherence and non-linear atom-optics of condensates. The 
remarkably powerful techniques for manipulating atomic condensates perfected in recent 
years suggest experiments difficult to imagine in the context of He superfluids or supercon- 
ductors, such as those now common-place in nonlinear and quantum optics, from four-wave 
mixing, parametric amplification, squeezed states, to Quantum Electrodynamics of cavity 
Bose condensates. 
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